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The  data  from  two  high-frequency  radars  (HFR)  operating  in  the  Iroise  Sea  are  re-processed  by  applying  an 
improved  version  of  the  direction  finding  algorithm,  removing  wave-induced  surface  currents  and  the  vari¬ 
ational  interpolation  on  a  regular  grid.  Combining  these  processing  techniques  allowed  reconstruction  of 
the  surface  currents  at  a  level  of  details  that  was  not  previously  available.  Refined  resolution  enabled  to  iden¬ 
tify  fine-scale  structures  of  surface  circulation,  to  quantify  the  variability  of  tidal  currents  and  the  residual 
(time  averaged)  velocity  field,  and  to  explain  spatial  intermittence  in  polarization  of  the  tidal  current  ellipses. 
The  analyzed  data  span  two  month-long  periods  in  spring  and  late  summer  of  2007.  The  major  findings  in¬ 
clude  (a)  a  dipole  structure  in  the  vorticity  field  characterized  by  two  oppositely  rotating  eddies,  generated 
on  the  leeward  side  of  the  Ushant  Island  at  flood  (negative  polarity)  and  at  ebb  (positive  polarity);  (b)  an  ex¬ 
tremely  strong  fortnightly  variability  of  tidal  currents  northwest  of  the  Ushant  Island  with  the  highest  veloc¬ 
ity  magnitude  of  3.9  m/s  caused  by  the  interference  of  the  major  semi-diurnal  tidal  constituents;  (c)  a 
significant  contribution  of  the  higher  order  nonlinear  tidal  harmonics  to  the  surface  currents  in  the  Fromveur 
strait,  which  maintains  strong  tidal  currents  and  affects  the  shape  of  their  fortnightly  modulation.  The  resid¬ 
ual  circulation  is  characterized  by  two  distinct  zones  approximately  separated  by  the  100  m  isobath:  in  the 
offshore  zone  the  residual  currents  have  a  significant  contribution  of  the  wind-driven  component,  whereas 
the  nearshore  zone  is  characterized  by  extremely  strong  (up  to  0.4  m/s)  time-independent  residual  circula¬ 
tion  featuring  two  permanent  anticyclonic  eddies:  north  of  the  western  extremity  of  the  Sein  archipelago, 
and  north  the  Ushant  Island.  The  acquired  data  and  the  presented  results  could  be  useful  for  regional 
model  validation  and  studies  of  the  local  eddy  dynamics,  tidal  fronts,  and  passive  tracer  transport  in  the 
region. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  high-frequency  (HF)  Doppler  radar  systems  have 
had  stunning  success  in  the  mapping  of  surface  currents.  The  ability 
to  map  surface  circulation  in  coastal  ocean  areas  has  brought  new  in¬ 
sights  to  the  complexities  of  physical  processes  in  nearshore  waters, 
and  allowed  significant  advances  in  our  understanding  of  circulation 
and  oceanographic  conditions  in  many  coastal  regions  (e.g.  Bassin  et 
al.,  2005;  Breivik  and  Saetra,  2001;  Haus  et  al.,  2000;  Kaplan  et  al., 
2005;  Kovacevic  et  al.,  2004;  Marmorino  et  al.,  1999;  Sentchev  et  al., 
2009a,  2009b;  Yoshikawa  et  al.,  2007). 

Iroise  Sea  is  one  of  the  regions,  continuously  monitored  by  HF  ra¬ 
dars  since  2006  (Fig.  1).  A  system,  composed  of  two  high-frequency 
Wellen  Radars  (WERA)  operating  at  12.4  MHz,  is  deployed  on  the 
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western  Brittany  coast  to  monitor  surface  currents  up  to  140  km  off¬ 
shore.  The  Iroise  Sea  circulation  is  very  difficult  to  study  by  in  situ 
methods,  as  it  is  dominated  by  extremely  strong  tidal  currents, 
often  exceeding  3  m/s,  and  affected  by  strong  and  often  violent  west¬ 
erly  and  southwesterly  winds.  Low  pressure  atmospheric  systems 
(cyclones),  generated  in  the  Northwestern  Atlantic,  regularly  cross 
the  Iroise  Sea  and  might  cause  significant  surges  on  the  West  Brittany 
coast  with  sea  surface  rise  up  to  0.50  m  (Bouligand  and  Pirazzoli, 
1999).  Surface  currents  in  the  Iroise  Sea  are  also  affected  by  swell 
and  wind  waves.  Violent  storms,  often  occurring  in  winter,  contribute 
to  the  reputation  of  the  Iroise  Sea  as  one  of  the  most  dangerous  seas 
in  Europe. 

Circulation  in  the  Iroise  Sea  has  been  studied  previously  by  many 
authors  by  means  of  field  measurements  (Le  Boyer  et  al.,  2009; 
Mariette  and  Le  Cann,  1985;  Pingree  et  al.,  1975)  and  numerical 
modelling  (Mariette  et  al.,  1982;  Muller  et  al.,  2009,  2010).  These 
studies  investigated  the  character  of  local  geostrophic  and  tidal  circu¬ 
lations,  eddy  frontal  dynamics  and  shed  light  on  the  role  of  tidal  mix¬ 
ing  in  controlling  cross-frontal  heat  and  mass  exchange.  They  also 
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Distance  (km) 


Fig.  1.  Location  of  the  Iroise  Sea  and  HFR  survey  area  (grey  square  contour  on  the  right  panel).  Experimental  domain  in  the  Iroise  Sea  (grey  shading)  and  radar  coverage  zone  (left 
panel).  Radar  sites  are  shown  by  grey  circles.  Grey  broken  contours  show  coverage  area  for  individual  radars.  Squares  denote  locations  of  the  met  buoy  and  tidal  gauge  station.  Open 
circle  indicates  the  location  of  the  grid  point  where  wave  induced  component  of  surface  current  was  analyzed  in  detail  (cf.  Fig.  2  d  and  Fig.  3b).  Contour  interval  of  the  bathymetry  is 
50  m  (grey  solid  lines).  Grey  dashed  lines  show  coverage  areas  for  the  individual  radars.  Geographic  names  used  in  the  text  are  also  shown. 


allowed  better  understanding  of  the  role  that  different  physical  forc¬ 
ing  mechanisms  play  in  the  Lagrangian  circulation  of  the  basin.  In 
particular,  the  influence  of  the  bottom  thermal  front  on  surface  circu¬ 
lation  and  large-scale  eddy  motions  associated  with  frontal  dynamics 
were  evidenced  in  (Le  Boyer  et  al.,  2009;  Muller  et  al.,  2009).  Howev¬ 
er,  with  the  exception  of  the  recent  field  study  of  Le  Boyer  et  al. 
(2009),  it  is  difficult  to  find  a  detailed  observation-based  description 
of  the  circulation  in  this  area  in  literature.  Quantitative  estimates  of 
the  small  scale  circulation  features,  their  magnitude  and  dependence 
on  various  forcing  factors  remain  relatively  uncertain.  This  shortcom¬ 
ing  is  primarily  due  to  the  difficulties  in  acquisition  of  the  long-term 
high-resolution  in  situ  data  because  of  the  extremely  strong  currents 
and  severe  meteorological  conditions  in  this  part  of  the  Atlantic 
Ocean.  Complexity  of  the  physical  processes  occurring  in  the  Iroise 
Sea  in  conjunction  with  complex  topography  renders  modelling  of 
the  circulation  also  difficult.  In  that  respect,  remote  sensing  of  surface 
currents  by  HF  radars  (HFR)  provides  a  unique  opportunity  to  estab¬ 
lish  a  monitoring  system  in  the  basin  on  a  regular  basis. 

As  the  experimental  data  required  for  model  validation  have  been 
extremely  scanty  until  recently,  in  2005,  the  French  Naval  Oceano¬ 
graphic  Centre  (SHOM)  put  forward  an  initiative  to  perform  the 
HFR  observations  of  surface  circulation  in  the  Iroise  Sea.  The  experi¬ 
ment  was  designed  to  get  better  knowledge  on  physical  processes 
governing  the  circulation  and  to  perform  the  validation  of  numerical 
models.  Beam  Forming  (BF)  algorithm  actually  used  for  processing 
the  radar  data  in  nearly  real-time  (Mariette  et  al.,  2006)  provides  ra¬ 
dial  velocities  of  surface  current  along  sixteen  beams  with  9°  azi¬ 
muthal  discretization.  The  angular  resolution  of  the  system  is  too 
coarse  and,  to  our  opinion,  unable  to  capture  the  fine  scale  structure 
of  the  surface  flow,  especially  around  the  islands  and  in  the  Fromveur 
Strait  which  are  of  particular  interest. 

In  the  present  study  we  address  two  challenges. 


First,  we  perform  high  resolution  surface  current  mapping  in  both 
temporal  and  spatial  dimensions  (20  min  and  1  km)  by  applying  a 
novel  approach  to  the  raw  HFR  data  processing  in  conjunction  with 
high-performance  variational  interpolation  technique.  High  azimuth¬ 
al  resolution  of  the  radial  velocities  is  achieved  by  the  HFR  signal  pro¬ 
cessing  using  an  improved  version  of  the  MUSIC  direction  finding 
algorithm  (Lipa  et  al.,  2006).  The  radial  velocities  are  then  interpolat¬ 
ed  on  the  regular  grid  using  a  variational  2dVar  method  (Yaremchuk 
and  Sentchev,  2009,  referred  hereafter  as  YS09).  Furthermore,  the 
contribution  of  surface  waves  to  the  radar  derived  velocities  is  esti¬ 
mated  with  the  method  recently  developed  by  Ardhuin  et  al. 
(2009).  These  three  techniques  have  been  implemented  for  proces¬ 
sing  of  HFR  data  for  two  selected  periods:  in  spring  and  late  summer 
2007  representing  in  total  52  days.  The  summer  period  is  synchro¬ 
nized  with  the  extensive  field  experiment  in  the  Iroise  Sea  abroad  of 
the  R/V  “Cotes  de  la  Manche”  (Le  Boyer  et  al.,  2009). 

Second,  we  study  the  resulting  surface  circulation  patterns  from 
the  perspective  of  different  modes  of  variability.  We  focus  on  physical 
processes  governing  the  tidal  currents  and  residual  circulation  in  the 
area.  The  HFR  observations  revealed  flow  complexity  and  variability 
at  a  level  of  details  that  were  not  previously  available.  Refined  resolu¬ 
tion  allowed  us  to  identify  fine-scale  structures  of  surface  circulation, 
to  quantify  the  variability  of  tidal  currents  and  residual  velocity  field, 
and  to  explain  patchiness  in  polarization  of  tidal  current  ellipses. 

The  paper  is  organized  as  follows.  In  Section  2,  we  present  the 
study  site,  the  environmental  data  (recorded  and  derived  from  nu¬ 
merical  modeling)  during  selected  periods  of  the  radar  experiment, 
the  HFR  network,  and  briefly  describe  the  methods  of  analysis.  In 
Section  3,  we  describe  and  analyze  the  observed  features  of  surface 
circulation,  flow  variability,  and  the  fine-scale  structures  in  vorticity 
and  residual  velocity  fields  in  the  Iroise  Sea.  Discussion  of  the  results 
and  summary  are  presented  in  Section  4. 
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2.  Data  and  methods 

2.1.  Study  site  and  environmental  conditions 

The  study  site,  Iroise  Sea,  is  located  in  the  extreme  north-eastern 
part  of  the  Gulf  of  Biscay  (Fig.  1).  This  is  a  shallow  water  area  with 
depth  gradually  decreasing  westward  from  50  to  150  m  with  more 
gentle  bottom  slope  in  the  South.  In  the  North,  there  is  a  group  of 
small  islands,  islets  and  rocks  that  form  Molene  archipelago,  and  a 
bigger,  8x4 km,  Ushant  Island.  It  is  separated  from  the  Molene 
Islands  by  the  2  km  wide  and  50  m  deep  Fromveur  strait.  In  the 


southern  part,  there  is  a  group  of  islets  and  rocks  belonging  to  the 
Sein  archipelago.  Only  a  part  of  it  is  exposed  at  low  tide. 

In  the  Iroise  Sea,  tidal  currents  are  generally  characterized  by 
clockwise  rotation  in  response  to  the  joint  effect  of  the  Coriolis  force 
and  sea  surface  gradients  related  to  tidal  wave  propagation.  Incident 
tidal  waves  arrive  from  the  Gulf  of  Biscay,  travel  around  the  W.  Brit¬ 
tany  Peninsula  and  enter  the  English  Channel.  The  tidal  waves’  inter¬ 
actions  is  the  dominant  factor  that  determines  variability  of  the  sea 
surface  height  (SSH)  and  currents  in  the  region.  The  SSH  variability 
in  the  Le  Conquet  harbor  exposes  variations  of  2  to  7  m  (Fig.  2a), 
with  the  predominant  semi-diurnal  period,  small  diurnal  inequality 
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Fig.  2.  External  forcing  data,  (a)  Tidal  elevation  observed  in  Le  Conquet  harbour  in  April-May  2007.  Wind  measurements  made  at  Met  Buoy  within  the  radar  coverage  zone  during 
April-May  (b)  and  August-September  2007  (c)  (Regional  Met-office  data).  Arrows  show  observed  wind  direction.  Wind  speed  scale  is  given  in  the  right  upper  corner,  (d)  3-h  data 
of  significant  wave  height,  Hs  (solid  line),  and  wave-induced  velocity,  Uss  (dashed  line),  in  August-September  2007  in  the  middle  of  the  study  area.  See  Fig.  1  for  location  of  mea¬ 
surement  points. 
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and  pronounced  fortnightly  modulation  due  to  the  interference  of  the 
semi-diurnal  (M2,  S2,  N2)  and  quarter-diurnal  (M4,  MS4)  constituents. 
Both  primary  and  secondary  spring  tides  occur  during  the  period 
shown.  Tidal  velocities  and  transports  show  fortnightly  variability  in 
response  to  the  spring-neap  cycle. 

Buoyancy  forcing  in  the  domain  is  relatively  weak.  The  majority  of 
freshwater  comes  from  the  Aulne  and  smaller  rivers  located  on  the 
western  Brittany  coast.  Their  discharge  to  the  Iroise  Sea  of  approxi¬ 
mately  60  m3/s  observed  in  April  2007  (http://www.hydro. 
eaufrance.fr)  can  be  considered  as  negligible,  compared  to  discharge 
of  the  Loire  and  Garonne  rivers  into  the  central  part  of  the  Gulf  of 
Biscay  (1200  and  900m3/s  respectively)  during  the  same  period. 
Great  distance  between  these  rivers'  mouths  and  the  Iroise  Sea  favors 
the  fresh  water  dispersion  and  only  a  small  quantity  of  freshwater 
reaches  the  western  Brittany  coast. 

Fig.  2b  shows  winds  measured  at  the  offshore  met  buoy  (see  Fig.  1 
for  buoy  location).  Two  distinct  wind  regimes  were  observed  during 
the  spring  period.  Between  the  10th  and  22th  of  April,  (the  first  12- 
day  period),  winds  were  generally  coming  from  the  North  and 
North-East  with  moderate  speeds  ranging  from  5  to  8  m/s.  From 
April  22  until  May  5,  variable  wind  regime  was  observed  with  small 
domination  of  northeasterly  winds  with  speeds  up  to  10  m/s.  During 
the  late  summer  period,  the  winds,  generally  stronger,  were  blowing 
mostly  from  the  northern  and  north-eastern  sectors  (Fig.  2c). 

The  spectrum  of  the  hourly  winds  recorded  by  the  offshore  Met 
buoy  is  given  in  Fig.  3a.  Spectral  analysis  of  both  zonal  and  meridional 
wind  components  and  the  total  kinetic  energy  of  the  wind  shows  the 
red  nature  of  the  spectrum  which  follows  a  power  low  with  spectral 
slope  close  to  —5/3,  mainly  in  higher  frequency  band  (>1  cyc/d). 
The  spectrum  reveals  a  variety  of  atmospheric  motions  with  signifi¬ 
cant  peaks  at  diurnal  and  sub-diurnal  frequencies  (0.55,  0.35,  and 
0.2  cyc/d).  Diurnal  peak  evidences  the  effect  of  sea  breeze,  whereas 
other  significant  peaks  are  within  the  range  of  cyclonic-anticyclonic 
frequencies.  They  are  certainly  related  with  cyclonic  activity  over 
the  North  Atlantic.  Small  scale  variability  of  the  local  winds  has 


been  documented  by  Muller  et  al.  (2007)  who  found  a  slight  system¬ 
atic  variation  of  the  wind  direction  in  the  area  between  the  Ushant 
Island  and  the  Brittany  coast. 


2.2.  HFR  data  and  methods  of  processing 
22  A.  HFR  velocity  data 

A  system  of  two  high-frequency  Wellen  Radars  ( WERA)  operating 
at  12.4  MHz  is  deployed  on  the  western  Brittany  coast  from  July  2006 
to  date.  Individual  radar  sites  are  located  at  Cape  Garchine  (site  G),  a 
seashore  flat-ground  area,  and  Cape  Brezellec  (site  B),  50  km  south¬ 
ward.  At  both  sites,  transmission  is  done  by  an  endfire  array  of  four 
antennas  forming  a  rectangle.  The  receiving  array  contains  16  equally 
spaced  antennas  parallel  to  the  shore  line  and  stretched  along  150  m. 

The  Beam  forming  (BF)  method,  actually  used  for  processing  the 
radar  data  in  near  real-time,  provides  radial  velocities  of  surface  cur¬ 
rent  along  beams  with  3  dB  width  of  9°  on  the  average.  An  automatic 
RFI  (Radio  Frequency  Interference)  clearing  has  been  integrated  into 
signal  processing  algorithm  (Gurgel  and  Barbin,  2008),  along  with 
the  quality  control  for  removing  outliers  and  radial  velocities  exceed¬ 
ing  4  m/s.  The  radial  (along  beam)  resolution  is  1.5  km  and  the  range 
is  of  the  order  of  140  km.  In  practice,  for  optimal  implementation  of 
the  RFI  clearing  method,  the  data  were  processed  up  to  a  range  of 
120  km. 

Radial  velocities  from  two  radars  are  combined  to  provide  surface 
current  maps  at  20-min  acquisition  rate  (http://www.previmer.org/ 
observations/courants/radar_hf_iroise).  The  accuracy  of  the  radar- 
derived  velocities  has  been  estimated  by  SHOM  (Oceanographic  Divi¬ 
sion  of  the  French  Navy)  through  a  comparison  with  surface  drifters 
and  ADCP  current  measurements  for  a  period  of  7  months.  In  the  ma¬ 
jority  of  situations,  the  discrepancy  in  velocity  measured  by  different 
instruments  did  not  exceed  0.15  m/s  (Le  Boyer  et  al.,  2009).  It  was 
also  noticed  that  a  large  fraction  of  this  discrepancy  might  not  be 
due  to  instrumental  errors,  but  to  the  definition  of  the  “surface 


Frequency  (cycles/day)  Frequency  (cycles/day) 

Fig.  3.  Power  spectra  of  the  wind  and  HF  radar  data,  (a)  Wind  spectrum  recorded  at  the  Met  buoy  in  April-May  2007.  (b)  Kinetic  energy  spectra  at  the  same  period  for  two  selected 
grid  points  of  the  radar  coverage  zone:  a  point  closest  to  the  Met  bouy  (black  line)  and  a  point  in  the  central  part  of  the  zone  (grey  line).  See  Fig.  1  for  point  location.  Numbers  stand 
for  dominant  tidal  period  (in  days).  Error  bar  indicates  95%  confidence  iterval. 
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velocity”,  as  other  measurement  techniques  are  sensitive  to  the  shear 
of  the  Eulerian  current,  to  Stokes  drift,  to  the  sea  state,  etc. 

Despite  a  good  quality  of  the  radial  velocities  recorded  by  the  ra¬ 
dars,  angular  resolution  of  the  system  is  relatively  coarse  and  does 
not  capture  the  fine  scale  structure  of  the  surface  flow,  especially 
around  Ushant-Molene  Islands  and  in  the  Fromveur  strait  which  are 
of  particular  interest. 

2.2.2.  Application  of  the  direction  finding  method  for  processing  the 
radar  data 

To  obtain  higher  resolution  radial  velocity  maps,  the  HFR  signals 
were  reprocessed  by  using  the  Multiple  Signal  Classification 
(MUSIC)  direction  finding  (DF)  algorithm  of  Schmidt  (1986)  driven 
by  the  parameterization  of  Lipa  et  al.  (2006).  In  contrast  to  BF,  the 
DF  method  seeks  for  the  bearings  from  which  the  energetic  Doppler 
lines  originate.  If  no  solution  is  delivered  in  a  given  direction  slot, 
then  the  DF  method  leaves  a  gap  in  the  velocity  record.  If  multiple 
Doppler  shift  solutions  are  delivered  in  the  same  bin,  a  power  aver¬ 
aged  Doppler  shift  value  is  calculated.  The  maximum  value  of  radial 
velocity  measured  by  a  radar  operating  at  a  wavelength  X 
(12.4  MHz)  is  limited  by/B-\/2  (4.3  m/s),  where  fB  is  the  Bragg  fre¬ 
quency  (0.359  Hz)  corresponding  to  the  radar  frequency. 

Processed  directions  spanned  the  aperture  from  —  70°  to  70°  with 
respect  to  the  normal  line  of  the  receiving  antenna  array,  thus  provid¬ 
ing  the  azimuthal  resolution  of  2°.  The  azimuthal  resolution  varies 
along  a  given  beam  from  0.35  km  (1.6  km  for  BF)  to  2.8  km 
(12.6  km  for  BF)  at  10  and  80  km  range,  respectively. 

Fig.  4  shows  radial  velocities  provided  by  BF  and  DF  methods  for 
the  same  acquisition  time.  In  general,  similar  spatial  patterns  on 
both  velocity  maps  are  found.  However  the  pattern  on  the  BF  map 
is  smeared  in  azimuth  and  it  can  be  difficult  to  determine  the  spatial 
extension  or  the  limits  for  different  structures  in  the  velocity  field. 
Detailed  analysis  of  velocity  maps  for  different  periods  of  tidal  cycle 
revealed  that  the  magnitude  and  spatial  variation  of  the  velocity 
field  were  poorly  represented  by  BF  algorithm.  However,  the  en¬ 
hanced  resolution  provided  by  the  DF  method  inherently  introduces 
gaps  (see  above)  and  noise  in  the  radial  velocity  field  (Fig.  4b). 
Noise  comes  from  the  fact  that  the  velocity  estimated  at  a  given 
range  for  a  particular  bearing  can  originate  from  another  bearing 
depending  on  the  complexity  of  the  ocean  wave  and  current  condi¬ 
tions  (Barrick  and  Lipa,  1996,  1997).  The  variational  method  used 
for  radial  velocity  interpolation  from  both  radar  units  (Section  2.2.3) 
performs  gap  filling  and  offers  an  efficient  control  of  the  noise  level. 


The  DF  technique  was  applied  to  the  HFR  data  to  retrieve  radial 
velocities  for  two  selected  periods,  namely,  from  10  April  to  5  May, 
and  from  24  August  to  19  September  2007.  The  second  period  was 
of  particular  interest  because  it  coincided  with  extensive  field  ex¬ 
periment  (Le  Boyer  et  al.,  2009).  During  these  periods,  both  radars 
worked  continuously  with  low  data  loss.  Fractional  availability  of 
radial  velocity  data  for  both  radars  during  the  two  periods  is 
shown  in  Fig.  5.  The  data  were  available  in  at  least  50%  of  time 
within  a  cone  with  a  radius  of  approximately  100  km.  Despite  of 
its  low  altitude  (20  m),  the  Ushant  Island  strongly  affects  availabil¬ 
ity  of  radar  data  providing  drops  in  the  data  return  to  25%  (for  site 
B)  and  to  50%  (for  site  G)  behind  the  island.  The  fractions  of  the 
data  return  by  the  individual  radars  are  shown  in  Fig.  5  by  thin  con¬ 
tours.  The  interpolation  grid  given  by  the  bold  contour  has  approx¬ 
imately  4000  points,  equally  spaced  by  1.1  km,  and  covers  a  nearly 
circular  area  with  the  seaward  limit  of  90  km  from  the  coast.  The 
condition  of  50%  availability  of  the  radial  velocity  data  is  met  in 
the  majority  of  the  interpolation  area  with  the  exception  of  the 
northern  part  for  site  B.  The  angle  between  radar  beams  at  points 
of  beam  intersection  was  greater  than  30°  all  over  the  interpolation 
area. 


2.2.3.  Variational  interpolation 

Velocity  vector  maps  were  generated  from  radial  velocity  data 
using  the  variational  interpolation  technique,  2dVar  (YS09).  This  is 
non-local,  kinematically  constrained  algorithm.  It  uses  a  combination 
of  radial  velocities  at  all  measurement  points  to  reconstruct  the  veloc¬ 
ity  vector  in  one  location.  This  feature  of  the  interpolating  technique 
can  help  to  overcome  some  limitations  related  to  a  lack  of  data.  The 
smoothness  of  the  interpolated  vector  field  is  enforced  by  penalizing 
the  squared  Laplacian  of  velocity  and  also  of  its  curl  and  divergence, 
thus  providing  the  algorithm  with  flexibility  and  efficient  control  of 
smoothness.  Kinematic  constraints  (zero  flux  on  rigid  boundary),  in¬ 
corporated  into  2dVar,  appear  particularly  useful  given  a  large  num¬ 
ber  of  islands  in  the  study  area.  The  2dVar  interpolation  was 
developed  and  employed  to  reconstruct  vector  current  fields  using 
synthetic  data  with  different  (up  to  50%)  signal/noise  ratios  and  also 
a  real  data  set  from  the  radar  experiment  in  Bodega  Bay,  northern  Cal¬ 
ifornia,  in  the  Gulf  of  Lion,  western  Mediterranean  (Sentchev  et  al., 
2009a,  2009b),  and  demonstrated  certain  advantages  compared  to  al¬ 
ternative  interpolation  techniques  (e.g.,  local  interpolation,  OMA 
(Muller  et  al.,  2009)). 
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Fig.  4.  Radial  velocity  maps  on  April  17,  2007  (4:10  am)  derived  from  radar  data  processing  by  (a)  beam  forming  and  (b)  direction  finding  methods. 
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Fig.  5.  Cumulative  availability  of  the  radial  velocity  data  (in  per  cent)  for  radar  station  B  (a)  and  G  (b)  during  the  two  periods  of  analysis  (thin  contour  lines).  Bold  contour  envelops 
the  area  where  surface  velocity  vectors  were  generated  at  each  time  interval  by  the  variational  algorithm. 


2.2.4.  Velocity  errors  estimation 

The  variational  approach  used  for  velocity  interpolation  allows  es¬ 
timating  errors  of  the  velocity  field  (YS09).  The  inverse  error  covari¬ 
ance  matrix  C-1  for  gridded  velocity  is  the  matrix  H  of  second 
derivatives  of  the  cost  function  with  respect  to  u  and  v  (the  Hessian 
matrix  of  the  variational  problem).  As  the  matrix  H  may  not  necessar¬ 
ily  have  the  full  rank,  we  seek  for  its  generalized  inverse  by  discarding 
eigenvectors  corresponding  to  small  eigenvalues.  Typically,  0.02%  of 
the  eigenvectors  are  rejected  during  the  inversion  of  H. 

Fig.  6  shows,  as  an  example,  the  map  of  errors  of  interpolated  ve¬ 
locity  (current  vector  magnitude)  for  the  first  HFR  snap  shot  (00:10 
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Fig.  6.  Velocity  errors  (cm/s)  for  a  surface  current  vector  map  on  April  10,  2007,  00:10. 


on  April  10,  2007).  The  map  combines  information  about  the  avail¬ 
ability  of  the  data,  accuracy  of  radial  velocity  measurements,  and  geo¬ 
metric  dilution  of  precision  (GDOP).  The  resultant  velocity  errors  are 
less  than  15  cm/s  in  the  majority  of  the  domain,  and  rapidly  increase 
(up  to  40  cm/s)  in  areas  poorly  covered  by  observations.  The  condi¬ 
tion  on  angle  intersection  (<30° ),  adopted  for  configuration  of  the  in¬ 
terpolation  domain,  prevents  error  growth  at  far  ranges.  It  should  be 
emphasized  that  the  error  distribution  corresponding  to  velocities  of 
the  observed  circulation  pattern  does  not  give  the  observational  sig¬ 
nal  to  noise  ratio.  This  ratio  can  be  derived  by  comparing  the  radial 
velocity  maps  with  spatial  distribution  of  the  accuracy  of  the  radial 
velocity  measurements.  The  latter  depends  upon  the  system  configu¬ 
ration,  level  of  backscattered  energy,  sea  surface  state,  variability  of 
currents  within  the  radar  measurement  cell,  ...,  and  might  attain 
8-10  cm/s. 

2.2.5.  Removing  the  sea  state  influence 

Surface  velocity  derived  from  radar  data  is  estimated  from  the  dif¬ 
ference  between  the  position  of  the  observed  and  theoretical  Bragg 
lines  on  the  Doppler  spectrum  of  the  sea  surface  echo  (Stewart  and 
Joy,  1974).  It  is  generally  assumed  that  Bragg  waves,  i.e.  waves 
whose  wavelength  (at  grazing  angle)  is  equal  to  the  half  of  the  elec¬ 
tromagnetic  wavelength,  satisfy  the  linear  dispersion  relationship  of 
ocean  waves.  Considering  nonlinearities  of  the  ocean  surface,  Weber 
and  Barrick  (1977)  slightly  modified  this  relationship  accounting  for 
all  types  of  the  surface  gravity  waves.  Broche  et  al.  (1983)  provided 
a  theoretical  expression  which  quantifies  the  effect  and  can  be  used 
to  extract  the  wave  induced  current  component  from  current  radio 
measurements.  The  method  of  velocity  correction  was  recently  ap¬ 
plied  to  radar  measurements  in  the  Iroise  Sea,  with  a  particular  inter¬ 
est  of  studying  the  wind-current  relationship  (Ardhuin  et  al.,  2009) 
locally.  In  the  present  paper  we  apply  the  method  to  the  whole  area 
of  radar  coverage.  Information  on  surface  waves  was  provided  every 
3  h  by  the  numerical  wave  model  WAVEWATCH  III  at  1/30°  resolution 
(2.5  km)  using  the  implementation  described  in  Ardhuin  et  al. 
(2009).  Appendix  A  contains  the  description  of  the  method  applied 
for  estimation  of  the  wave  induced  surface  current  velocity  vector, 
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Uss,  at  the  grid  points  of  the  wave  model.  Resulting  values  of  Uss  were 
spatially  interpolated  on  the  regular  radar  grid,  then  interpolated  in 
time,  and  finally  subtracted  from  the  radar  derived  current  velocities. 
For  convenience,  Uss  will  be  referred  hereafter  as  Stokes  velocity  al¬ 
though  in  practice  it  represents  approximately  85%  of  the  model- 
estimated  Stokes  drift  velocity  (cf.  Appendix). 

Fig.  2  d  provides  estimated  Stokes  velocity  and  significant  wave 
height  during  the  second  (late  summer)  period  of  analysis  in  a 
point  located  in  the  middle  of  the  study  area  (Fig.  1).  Strong  correla¬ 
tion  is  observed  between  the  wind  speed  and  Uss.  Northern  winds 
with  high  speed,  observed  on  August  29,  September  3  and  18 
(Fig.  2c),  produce  surface  currents  with  magnitudes  greater  than 
0.10  m/s.  Winds  from  North-East  sector  tend  to  weaken  surface 
waves  (Hs<  1  m)  and  the  associated  wave  induced  currents  to  less 
than  0.10  m/s  (Fig.  2  d). 

Fig.  7  summarizes  the  contribution  of  Stokes  velocity  to  the  sur¬ 
face  currents  in  terms  of  the  occurrence  probability  of  a  measured 
current  magnitude  U  and  the  wave  induced  current  velocity  Uss . 
These  statistical  estimates  cover  both  studied  periods.  The  modulus 
of  Uss  is  a  small  quantity  compared  to  U  with  a  peak  of  occurrence 
at  0.05  m/s  whereas  the  occurrence  peak  for  U  is  0.40  m/s.  The  high¬ 
est  magnitude  of  Stokes  current  did  not  exceed  0.18  m/s  during  the 
period  of  interest.  Moreover,  Uss  values  are  rather  small  compared 
to  the  precision  of  radar  surface  currents.  For  a  point  located  in  the 
middle  of  the  study  area  (Fig.  1),  the  absolute  values  of  the  west- 
east  (south-north)  components  of  Uss  exceed  the  intrinsic  error  of 
the  corresponding  current  components,  estimated  as  4  cm/s  (7  cm/ 
s)  in  53%  (7%)  of  cases.  However  the  west-east  component  of  Uss  is 
generally  higher  than  this  error  and  can  exceed  it  4  times.  In  that 
case,  removing  the  sea  state  influence  is  useful  to  improve  the  quality 
of  current  velocity  estimates. 

2.2.6.  Principal  component  and  rotary  spectral  analysis 

In  order  to  identify  processes  governing  surface  circulation  in  the 
Iroise  Sea,  spectral  analysis  of  velocity  records  in  different  locations 
has  been  performed.  Fig.  3b  shows  spectral  distribution  of  the  total  ki¬ 
netic  energy  for  two  selected  grid  points:  a  point  closest  to  the  Met 
bouy  and  a  point  in  the  central  part  of  the  radar  coverage  zone  (see 
Fig.  1  for  locations).  The  distribution  clearly  shows  that  surface  dy¬ 
namics  in  the  Iroise  Sea  is  dominated  by  tides.  The  spectra  in  Fig.  3b 
indicate  that  most  of  the  kinetic  energy  is  distributed  between  the 
dominant  tidal  frequencies:  semi-diurnal,  quarter-diurnal,  sixth- 
diurnal  and  diurnal,  with  prevailing  semi-diurnal  tide.  In  sub-tidal 


velocity  m/s 

Fig.  7.  Probability  of  occurrence  of  the  velocity  magnitude  of  surface  currents  derived 
from  radar  measurements  and  the  velocity  magnitude  of  wave  induced  currents 
(Stokes  current). 


frequency  range,  there  is  a  moderate  rise  of  energy  which  reveals  mo¬ 
tions  with  3.3-day  period  (Fig.  3b).  Analysis  of  similar  distribution  for 
the  wind  (Fig.  3a)  does  not  show  any  significant  atmospheric  forcing 
in  tidal  band,  suggesting  the  absence  of  modulation  of  tidal  currents 
by  wind. 

Furthermore,  we  applied  the  principal  component  analysis  (PCA) 
technique  to  this  data  set  in  order  to  quantify  tidal  flow  dynamics. 
Our  approach  has  a  certain  advantage  over  a  frequently  used  harmonic 
analysis  because  it  allows  quantifying  the  total  contribution  of  all  tidal 
constituents  to  observed  currents  and  assessing  time-space  variability 
of  the  currents.  As  tidal  currents  in  the  Iroise  Sea  are  rotational,  the 
current  velocity  vector  evolving  over  a  tidal  cycle  draws  an  ellipse.  Pa¬ 
rameters  of  synthesized  ellipses,  retrieved  from  the  PCA,  provide  ori¬ 
entation  and  magnitude  of  the  dominant  current.  The  anisotropy  of 
oscillatory  tidal  flow  (ellipse  eccentricity)  is  quantified  by  estimating 
the  eigenvalue  ratio  of  the  velocity  correlation  tensor.  The  PCA  tech¬ 
nique  does  not  give,  however,  information  on  the  sign  of  rotation  of 
velocity  vectors.  To  fill  this  gap,  we  performed  rotary  spectral  analysis 
(Emery  and  Thompson,  1997)  of  velocity  time  series,  referred  hereaf¬ 
ter  as  RSA.  The  technique  involves  the  decomposition  of  the  velocity 
vector  into  the  clockwise  (cw)  and  counter-clockwise  (ccw)  rotating 
circular  complex-valued  components.  First,  the  RSA  was  performed 
to  identify  the  dominant  frequency  in  current  velocity  variation.  Both 
cw  S_  (negative)  and  ccw  S+  (positive)  power  spectra  revealed  pro¬ 
nounced  peaks  at  the  semi-diurnal  frequency.  After  that,  the  rotary  co¬ 
efficient,  r=  (S+  —  S_)/(S+  +S_),  was  estimated  at  every  grid  point  at 
these  peak  values,  r  ranges  from  —  1  for  clockwise  motion  to  + 1  for 
counter-clockwise  motion  (r=0  is  oscillating  non-rotational  flow). 

We  have  also  estimated  the  Eulerian  residual  currents  and  investi¬ 
gated  their  relationship  with  wind  forcing. 


3.  Results 

3.1.  Tidal  currents:  spatial  patterns  and  variability 

Spatial  and  temporal  variability  of  tidal  currents  is  quantified 
using  the  parameters  of  synthesized  tidal  current  ellipses  derived 
from  the  PCA  and  RSA.  The  results  are  summarized  in  Fig.  8  for  the 
primary  spring  tide  period.  Tidal  current  ellipses  (Fig.  8a)  represent 
the  7-day  averaged  circulation  pattern  associated  with  the  tidal 
wave  arriving  from  the  Gulf  of  Biscay  and  traveling  northeastward  to¬ 
ward  the  English  Channel.  The  detailed  analysis  of  ellipse  orientation 
shows  that  tidal  currents  are  controlled  by  the  bottom  topography 
and  peculiarities  of  the  coastline.  These  factors  tend  to  align  the 
major  axes  along  the  depth  contours. 

One  can  distinguish  anisotropy  of  the  current  field  with  relatively 
high  ellipse  eccentricity  in  the  western  part  of  the  region  and  in  two 
locations  in  the  eastern  part:  south  of  the  Ushant  Is.  and  in  a  shallow 
water  region  between  the  Molene  Islands  and  Brittany  coast.  The  for¬ 
mer  might  be  related  to  the  current  squash  against  the  island  produc¬ 
ing  more  circular  shape  of  ellipses.  The  latter  reflects  the  effect  of 
bottom  friction  on  tidal  currents  (Defant,  1961). 

The  spatial  distribution  of  the  rotary  coefficient  (Fig.  8b)  shows 
that  currents  are  rotating  clockwise  (cw)  in  the  majority  of  the  do¬ 
main.  This  feature  is  related  to  the  spatial  pattern  of  the  sea  surface 
gradients  which  tend  to  follow  the  propagation  of  the  tidal  wave 
from  the  Gulf  of  Biscay  in  a  broad  cw  turn  around  the  Brittany  Penin¬ 
sula  (e.g.  Lux  et  al.,  2010).  The  Coriolis  force  acting  on  the  currents  re¬ 
inforces  their  cw  rotation  (Defant,  1961). 

The  maximum  negative  values  of  r  (up  to  —0.9)  are  found  in  the 
northeastern  part  of  the  region  (Fig.  8b).  Large  positive  values  (up 
to  0.7)  are  observed  in  three  locations  around  the  Ushant  and  Molene 
Islands  (Fig.  8b),  and  in  the  vicinity  of  the  Sein  archipelago  (not 
shown).  In  the  shallow  water  areas,  east  of  Molene  and  north  of 
Sein  islands,  the  bottom  friction  comes  into  play,  alternating  the 
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patchiness  in  rotary  coefficient  distribution  with  change  in  sign  of  ro¬ 
tation  of  the  velocity  vectors  around  the  islands  is  the  most  notable 
feature  of  tidal  circulation  in  the  Iroise  Sea  (Fig.  8a, b). 

The  magnitude  of  tidal  currents  shows  significant  variations  in 
space  and  time.  During  spring  tide,  current  velocity  varies  from  0.7 
to  3.9  m/s  (Fig.  8c),  during  neap  tide  —  from  0.2  to  1.7  m/s  (results 
not  shown). 

Tidal  waves  arriving  from  the  Gulf  of  Biscay  propagate  in  the  N-NE 
direction  over  steep  bottom  in  the  southern  part  of  the  Iroise  Sea,  and 
spring  tide  velocities  do  not  exceed  1.5  m/s  there.  Higher  bathymetry 
gradients  in  the  NE  part  of  the  Sea  and  the  presence  of  islands  cause 
flow  acceleration,  tend  to  tighten  the  streamlines  and  provide  the 
maximum  velocities  at  the  seaward  extremity  of  the  Ushant  Is.  (up 
to  3.9  m/s)  and  in  the  Fromveur  Strait  (up  to  3.8  m/s).  We  performed 
detailed  analysis  of  surface  current  variability  in  these  two  locations. 

Figs.  9a  and  10a  show  time  evolution  of  the  magnitude  of  tidal 
currents  (semi-major  axis  of  synthesized  current  ellipses)  along  two 
lines  crossing  the  region  with  strongest  currents.  The  first  line,  orient¬ 
ed  in  the  west-east  direction  (WE-line)  is  located  north  of  Ushant  Is. 
The  second  passes  through  the  Fromveur  Strait  (FS-line)  and  has  an 
oblique  orientation  (Fig.  8c). 

Both  primary  and  secondary  spring  tides  can  be  recognized  in  the 
velocity  field.  In  the  central  part  of  the  WE-line,  at  point  A  located  NW 
of  Ushant  Is.  (Fig.  8c),  the  magnitude  of  tidal  currents  increases  more 
than  four  times  in  8  days,  rising  from  less  than  1  m/s  on  1 1  April  to 
3.9  m/s  on  18  April  (Fig.  9b).  The  range  of  current  variability  de¬ 
creases  progressively  westward  and  eastward  from  the  maximum 
value  at  point  A.  The  ellipse  eccentricity  and  orientation  do  not 
change  with  time  (Fig.  9b). 

Variability  of  the  tidal  flow  in  the  Fromveur  Strait  (point  B)  looks 
quite  different  (Fig.  10b).  The  current  velocity  rises  from  1.7  to  3.8  m / 
s  during  the  same  period  whereas  the  curve  of  current  magnitude 
variation  with  time  in  this  point  and  in  the  neighboring  points  is  flat¬ 
tened  during  spring  tide,  suggesting  that  the  fortnightly  modulation 
of  the  velocity  is  affected  by  non-linear  interactions  between  the 
tidal  constituents.  To  explore  this  hypothesis,  we  performed  harmon¬ 
ic  analysis  of  tidal  currents  in  two  particular  locations:  north-west  of 
Ushant  Is.  (point  A),  and  in  the  Fromveur  strait  (point  B).  The  velocity 
time  series  were  analysed  using  the  Matlab  software  of  Pawlowicz  et 
al.  (2002).  The  results  reveal  noticeably  stronger  current  magnitude 
(semi-major  axis)  of  higher  order  tidal  harmonics  (M6,  2MS6, 
2SM6)  (Table  1)  indicating  their  amplification  in  the  strait.  This  fea¬ 
ture  also  reveals  complexity  of  tidal  dynamics  in  the  strait  and  dem¬ 
onstrates  that  shallow-water  constituents  might  play  an  important 
role  in  fortnightly  modulation  through  the  interaction  with  principal 
tidal  constituents. 


Fig.  8.  (a)  Tidal  current  ellipses  during  spring  tide  in  April  2007  (7-day  averaged  from 
15  to  22  April  2007).  Every  forth  ellipse  is  shown.  Red  ellipses  denote  counter¬ 
clockwise  rotating  currents,  blue  ones-clockwise  rotating  currents.  Rotary  coefficient 
r  for  the  same  period  (b),  and  amplitude  of  spring  tidal  currents  (c)  in  the  vicinity  of 
Ushant  and  Molene  Islands.  50  m  isobath  is  shown  by  thin  grey  contour.  Grey  lines  in 
(c)  show  cross-sections  used  for  analysis  of  the  variability  of  tidal  currents. 

sign  of  current  vector  rotation,  and  providing  positive  polarization  of 
current  ellipses. 

The  largest  region  of  the  ccw  rotating  currents  is  located  south  of 
Ushant  and  Molene  Islands.  Positive  polarization  cannot  be  caused  di¬ 
rectly  by  the  bottom  friction  as  the  depth  exceeds  50  m  there.  One  of 
the  possible  mechanisms  is  the  interaction  of  tidal  flow  with  the  Ush¬ 
ant  Is.  resulting  in  vorticity  production.  Positive  vorticity,  generated 
by  frictional  effects  on  leeward  side  of  the  island  during  the  south¬ 
ward  ebb  tidal  flow,  is  advected  westward  along  the  southern  periph¬ 
ery  of  the  Ushant  and  Molene  Islands  during  transition  from  ebb  to 
flood.  This  provides  an  extended  region  with  ccw  rotating  ellipses 
there  (Fig.  8a).  The  origin  of  the  negative  values  of  r  observed  north 
of  the  Ushant  Is.  (Fig.  8b)  can  be  explained  in  a  similar  way.  This 


3.2.  Vorticity  field 

Surface  current  vorticity,  averaged  over  the  25-day  period  in 
spring,  shows  patchiness  in  spatial  distribution  (Fig.lla)  and  some 
degree  of  similarity  with  the  spatial  pattern  of  rotary  coefficient 
(Fig.  8b).  While  in  the  majority  of  the  domain,  vorticity  does  not  ex¬ 
ceed  ±2  10-5s-1,  much  higher  values  are  found  around  the  islands 
and  in  the  south-eastern  sector.  These  extreme  values  reach  —18 
10-5s-1  for  negative,  and  15  10-5s-1  for  positive  vorticity 
(Fig.  11a)  thus  evidencing  a  pronounced  rotational  character  of  the 
surface  circulation  in  the  vicinity  of  islands. 

Analysis  of  instantaneous  velocity  patterns  allows  to  qualitatively 
explain  the  vorticity  distribution  and  to  assess  its  production.  The  se¬ 
lected  spring  tide  period  (18-19  April,  2007)  was  characterized  by 
moderate  northeasterly  winds  with  the  relevant  circulation  patterns 
shown  in  Fig.  llb-d.  The  ebb  spring  flow  is  affected  by  the  bottom 
friction  in  the  vicinity  of  islands.  Thus,  southward  tidal  currents  gen¬ 
erate  a  cyclonic  eddy  on  the  leeward  side  of  the  Ushant  Is.  (Fig.  lib) 
and  cause  positive  vorticity  production.  After  the  current  reversal 
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Fig.  9.  (a)  Time  evolution  of  the  amplitudes  of  tidal  currents  in  April-May  2007  along  the  WE-line  (Fig.  6c).  (b)  Evolution  of  tidal  current  ellipses  at  point  A  (Fig.  6c)  located  on  the 
WE-line  and  where  the  observed  surface  currents  are  highest. 


(Fig.  lie),  the  eddy  is  advected  westward  by  about  10-15  km  (loca¬ 
tion  of  the  second  positive  vorticity  maximum  in  Fig.  11a),  and  then 
disappears  during  the  developing  flood.  On  flood,  surface  currents  skirt 
around  the  Ushant  and  Molene  Islands  in  two  branches:  stronger  west¬ 
ern  stream,  and  near-shore  eastern  stream,  also  leaking  through  the 
Fromveur  Strait.  Frictional  effects  generate  negative  vorticity  north  of 
Ushant  and  Molene  Islands  and  form  an  eddy  (Fig.  lid).  This  ccw  rotat¬ 
ing  eddy  grows  until  the  current  reversal  at  high  water,  then  disappears 
during  the  development  of  the  ebb  flow. 

Formation  of  the  cw  eddy,  partially  captured  in  the  south-eastern 
sector  of  radar  coverage  zone  (results  not  shown)  follows  similar  sce¬ 
nario.  Negative  vorticity  is  generated  by  frictional  effects  of  Sein  ar¬ 
chipelago  on  the  oscillating  tidal  currents.  The  maximum  values  of 


vorticity  reach  —  1210_5s_1,  indicating  generation  of  an  anticyclon- 
ic  eddy  north  of  the  central  axis  of  the  archipelago. 

Another  specific  feature  of  vorticity  distribution  is  a  narrow  band  of 
negative  vorticity  surrounding  the  southwestern  and  eastern  (From¬ 
veur  St.)  parts  of  the  Ushant  Is.  (Fig.  11a).  This  vorticity  pattern  indicates 
the  dominant  cw  sign  of  surface  flow  around  the  islands.  A  transverse 
circulation  induced  by  a  tidal  current  escalating  the  depth  gradient 
was  analyzed  by  Godin  (1988)  who  showed  that  depth  gradients  and 
bottom  friction  tend  to  veer  the  mean  current  to  the  left  of  the  direction 
of  tidal  wave  propagation  (in  the  northern  hemisphere),  align  it  with 
the  depth  contours,  thus  producing  cw  rotation  of  the  mean  flow. 

Spatial  distribution  of  vorticity  in  summer  (results  not  shown) 
was  found  to  be  very  similar  to  that  observed  in  spring,  thus 
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Fig.  10.  (a)  Time  evolution  of  the  amplitudes  of  tidal  currents  in  April-May  2007  along  the  FS-line  (Fig.  6c).  (b)  Evolution  of  tidal  current  ellipses  at  point  B  (Fig.  6c)  located  in  the 
middle  of  the  Fromveur  strait  where  the  observed  surface  currents  are  highest. 


providing  additional  support  for  the  permanent  character  of  the  ob¬ 
served  circulation  (vorticity)  pattern  around  islands. 

3.3.  Residual  currents 

In  the  tidally  dominated  basins  of  the  Northwestern  European 
shelf  such  as  the  Iroise  Sea,  Ireland  Sea  or  the  English  Channel,  the 


Table  1 

Semi-major  axes  (m/s)  of  tidal  current  ellipses  for  the  most  important  tidal  constitu¬ 
ents  in  two  selected  locations:  north-west  of  the  Ushant  Island  (pt.  A)  and  in  the  From¬ 
veur  strait  (pt.  B  in  Fig.  6c). 


Tidal  constit. 

M2 

S2 

N2 

M4 

MS4 

2MN6 

M6 

2MS6 

A:  Ushant 

1.74 

0.92 

0.47 

0.15 

0.16 

0.06 

0.06 

0.08 

B:  Fromveur 

2.09 

0.83 

0.35 

0.12 

0.09 

0.12 

0.16 

0.20 

long-term  averaged  currents,  referred  hereafter  as  residual  currents 
(RC),  might  govern  the  transport  of  material  to  a  great  extent.  There¬ 
fore,  estimation  of  residual  circulation  is  important  for  assessing  the 
long-term  transport  and  dispersion  of  material  in  tidal  flows  (Bailly 
du  Bois  and  Dumas,  2005;  Foreman  et  al.,  1992;  Zimmerman,  1986). 
We  estimated  the  residual  velocities  by  averaging  the  monitored  sur¬ 
face  currents  over  each  of  the  two  periods  under  consideration.  The 
resultant  fields,  shown  in  Fig.  12a  and  b,  include  tidally  induced  com¬ 
ponent  (major  contribution),  the  wind  and  buoyancy  forced  residual 
flows.  The  wave  induced  (Stokes)  component  was  removed  before¬ 
hand  from  velocity  time  series.  Fig.  12  clearly  demonstrates  that  con¬ 
trol  of  the  residual  flow  by  tidal  motions  is  not  prevailing  over  the 
area  of  radar  survey,  and  that  tidal  and  wind  forcing  might  be  in 
competition. 

The  spatial  structure  of  residual  currents  for  both  periods  appears 
complex  and  allows  dividing  the  basin  into  two  zones;  the  nearshore 
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Fig.  11.  (a)  Average  vorticity  field  (x  10  5  s  *)  for  April-May  period  of  2007.  (b)  Surface  currents  around  Ushant  and  Molene  Islands  on  ebb  during  spring  tide.  Surface  currents  at 
the  beginning  (c)  and  at  the  end  of  flood  (d).  50  m  isobath  is  shown  by  thin  grey  contour. 


zone,  with  the  dominant  contribution  of  tides  to  RC,  and  the  offshore 
zone,  where  wind  effect  is  also  playing  a  significant  role.  These  two 
zones  can  be  separated  by  a  line  slightly  west  of  the  100-m  izobath. 
In  the  western  (offshore)  zone,  current  velocities  do  not  exceed 
0.15  m/s  in  April  and  0.25  m/s  in  September  with  westward  and 
south-westward  directions  respectively,  indicating  the  effect  of  dom¬ 
inant  northeasterly  and  northerly  winds  during  these  periods  (Fig.  2). 
On  the  contrary,  in  the  eastern  (nearshore)  zone,  the  RC  looks  similar 
for  both  periods  and  reveals  eddies  in  the  southeastern  sector  and 
north  of  the  Ushant  Is.  There  is  also  a  tendency  for  intensification  of 
the  nearshore  currents  (up  to  0.40  m/s)  along  the  50-m  izobath 
(Fig.  12).  The  second  region  of  elevated  velocities  is  aligned  along 
the  100-m  izobath.  Stronger  currents  (up  to  0.40  m/s)  occur  in  the 
south-eastern  part  of  the  domain  in  the  vicinity  of  Sein  archipelago. 
The  velocity  drops  to  0.20  m/s  further  northward,  and  the  flow 
seems  to  split  in  two  branches.  This  separation  and  deflection  of  cur¬ 
rents  westward  and  eastward  is  probably  caused  by  the  bottom  irreg¬ 
ularity  and  a  promontory  visible  on  the  bathymetry  map  at  the 
intersection  of  coordinate  lines  x  =  60  and  y  =  95  km  (Fig.  12a).  Cur¬ 
rents  entering  the  western  offshore  zone  change  direction  under  the 
effect  of  wind  forcing,  whereas  the  eastern  branch  merges  with  the 
coastal  current.  Although  the  residual  current  intensification  in  the 
Iroise  Sea  has  been  already  documented  (Mariette  et  al.,  1982; 
Muller  et  al.,  2010),  the  eddies  around  the  Ushant  Is.  in  the  RC  field, 
revealed  from  the  analysis  of  HFR  data  set,  is  a  novel  feature  of  tidal 
dynamics  in  the  area. 

The  maps  of  residual  velocities  for  both  periods  show  strong  wind 
forcing  control  in  the  western  zone  and  topographic  control  in  the 


eastern  zone.  The  response  of  surface  currents  to  wind  forcing  was 
further  examined  statistically  by  assuming  the  following  relation  for 
the  (non-tidal)  residual  times  series  (Prandle  and  Matthews,  1990): 

R(t)=R  +  ar(t-At),  (1) 

where  R  is  the  steady  non-wind-driven  residual  current,  At  is  a  time 
lag  between  wind  and  surface  current  response,  and  a  is  a  complex¬ 
valued  coefficient.  Both  R  and  a  vary  in  space,  r  is  the  wind  stress  cal¬ 
culated  from  the  hourly  wind  data  of  met  buoy  and  assumed  to  be 
power-law  function  of  the  wind  speed  W".  The  complex  current  vec¬ 
tor  R(t)  was  derived  from  low-pass  filtered  HFR  surface  current  veloc¬ 
ities.  The  filter  reduces  tidal  (diurnal,  semidiurnal)  and  inertial 
variability  to  less  than  1%  of  their  original  values,  and  passes  90%  or 
more  of  the  low-frequency  signal  (Stacey  et  al.,  1986).  The  wind 
hourly  data  were  smoothed  by  24-h  moving  average.  The  values  of 
R  and  a  were  calculated  by  least-squares  fitting  for  different  powers 
of  the  wind  speed,  W,  and  time  lag  At.  The  overall  correlations  be¬ 
tween  the  observed  and  estimated  residual  velocities  were  weakly 
sensitive  to  the  parameter  n  of  the  wind  speed,  but  showed  signifi¬ 
cant  dependence  on  the  time  lag  (in  the  range  0-12  h).  The  maxi¬ 
mum  correlation,  0.86,  was  found  for  a  time  lag  of  4  h  and  t  oc  W2. 
The  optimal  values  of  a  (providing  maximum  correlation)  were 
then  used  in  combination  with  the  wind  stress  to  subtract  the  wind 
response  from  surface  velocity  time  series. 

Two  situations  were  explored.  In  the  first  situation,  only  tidal  mo¬ 
tions  were  removed  from  velocity  time  series  R(t).  Then,  the  wave  in¬ 
duced  surface  currents,  retrieved  by  the  method  described  in 
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Fig.  12.  Eulerian  residual  surface  currents  during  the  26-day  period  in  April-May  2007 
(a)  and  late  summer  period,  August-September  2007,  (b).  Grey  lines  indicate  the  loca¬ 
tion  of  east-west  cross  section  used  for  assessment  of  residual  current  variability.  50  m 
and  1 00  m  isobaths  are  shown  by  thin  grey  contours. 


Section  2.2.5,  were  also  subtracted  from  these  time  series,  and  the 
least-squares  fit  of  surface  velocities  to  wind  forcing  was  performed. 
The  results  revealed  high  sensitivity  of  surface  residuals  to  wave- 
induced  forcing,  which  in  turn  contains  the  wind  effect.  The  contribu¬ 
tion  of  wave  induced  currents  to  surface  residuals  was  found  particu¬ 
larly  important  in  the  eastern  shallow  water  part  of  the  domain 
during  the  late  summer  period  (when  winds  were  blowing  from  the 
northern  sector).  The  values  of  wave-induced  velocity  in  the  vicinity 
of  Ushant  and  Molene  Islands  attain  0.15-0.18  m/s  (cf. 
Section  2.2.5)  and  become  comparable  with  the  magnitude  of  the  re¬ 
sidual  currents  in  this  area. 

In  Fig.  13a,  we  show  the  northern  component  of  de-tided  velocity 
time  series  for  April,  and  in  Fig.  13b,  velocities  from  which  the  major 
contribution  of  wind  and  waves  were  removed  by  linear  regression 
method.  We  examine  the  time  variability  of  this  velocity  component 
across  a  horizontal  line  located  approximately  in  the  middle  of  the 
HFR  coverage  zone  (Fig.  12).  Fig.  13a  reveals  existence  of  a  permanent 
northward  residual  current,  following  the  100-m  izobath.  The  veloci¬ 
ty  magnitude  varies  in  time  from  0.05  to  0.20  m/s  indicating  a  low- 
frequency  modulation  of  the  RC  by  atmospheric  forcing.  Northern 
and  north-eastern  winds  observed  during  the  first  4-day  period 
(Fig.  2b)  significantly  weaken  the  current  to  0.01  m/s.  In  the  eastern 
nearshore  part  of  the  section,  southward  residual  currents  persist 
during  the  whole  period  of  analysis,  with  velocities  up  to  0.10  m/s 
(Fig.  13a).  In  the  western  (offshore)  sector,  the  current  changes  the 
sign  to  positive  with  each  reversal  of  the  wind  direction. 

The  removal  of  wind-wave  effect  allows  better  understanding  of 
RC  variability  (Fig.  13b).  The  major  apparent  difference,  compared 


to  Fig.  13a,  is  the  dominant  northward  flow  over  the  majority  of  the 
cross-section,  with  the  exception  of  the  eastern  limit  influenced  by 
the  anti-cyclonic  eddy  (cf.  Fig.  12a).  The  change  in  sign  of  the  current 
direction,  observed  on  April  14-17  and  22-26,  indicates  a  low  fre¬ 
quency  variability  of  the  RC,  caused  by  a  mechanism  different  from 
wave  and  wind  forcing.  This  variability  also  results  in  a  periodic  in¬ 
crease  of  current  velocity  up  to  0.22  m/s  within  the  near-shore  jet. 

During  the  late  summer  period,  the  RC  looks  quite  different 
(Fig.  14):  the  northward  current  (approximately  following  the  100- 
m  izobath)  keeps  its  spatial  shape  and  high  magnitude  in  the  eastern 
part  of  the  section,  the  dominant  flow  appears  clearly  southward  in 
the  western  part  of  the  section.  Removal  of  the  wind-  and  wave- 
induced  currents  modifies  the  circulation  pattern  (Fig.  14b)  and  sug¬ 
gests  that  other  forcing  (e.g.  baroclinic)  might  come  into  play  and 
contribute  to  the  RC  formation  in  the  region.  Positive  velocities  of 
the  RC,  observed  close  to  the  western  limit  of  the  domain  between 
11  and  18  September  (Fig.  14b),  support  this  hypothesis:  field  studies 
performed  in  this  area  during  the  same  period  (Le  Boyer  et  al.,  2009), 
indicated  meandering  of  the  bottom  thermal  front  which  might  con¬ 
tribute  to  generation  of  a  filament  with  positive  (northward)  veloci¬ 
ties  separating  from  the  main  northward  branch  of  the  RC 
(Fig.  14b).  It  should  also  be  noted  that  a  signature  of  the  wind  forcing 
still  remains  non-negligible  in  the  RC  during  the  period  of  strong 
northern  winds.  Detailed  investigation  of  the  circulation  in  the  period 
of  Ushant  thermal  front  development  is  a  subject  of  our  future  work. 

4.  Discussion  and  conclusions 

The  experiment  reported  in  this  paper  is  based  on  the  employ¬ 
ment  of  a  HFR  system  to  study  circulation  in  highly  energetic  tidal 
basin  —  the  Iroise  Sea.  The  radar-derived  surface  currents  were  inves¬ 
tigated  for  two  selected  periods  by  means  of  various  numerical  and 
statistical  techniques.  We  managed  to  obtain  surface  currents  at 
higher  resolution,  due  to  application  of  the  MUSIC  (direction  finding) 
technique  in  conjunction  with  variational  2dVar  interpolation.  This 
approach  provided  significant  level  of  confidence  in  the  surface  flow 
retrieved  from  HFR  data.  Estimated  accuracy  of  the  surface  velocities 
was  found  to  be  less  than  0.15  m/s  (30%)  in  the  majority  of  the  radar 
coverage  zone  (Fig.  6),  whereas  the  mean  velocity,  averaged  in  space 
and  time  for  the  first  period  of  observations,  is  close  to  0.50  m/s. 

The  radar  derived  velocity  data  reveal  flow  complexity  and  vari¬ 
ability  at  a  level  of  details  that  were  not  previously  available.  Though 
there  have  been  a  number  of  important  field  investigations  (Ardhuin 
et  al.,  2009;  Le  Boyer  et  al.,  2009)  and  numerical  studies  of  circulation 
in  the  Iroise  Sea  (Muller  et  al.,  2009,  2010),  the  presented  radar  sur¬ 
vey  have  brought  a  significant  contribution  to  understating  physical 
processes  governing  the  coastal  water  dynamics  in  the  region  where 
the  tidal  currents  attain  one  of  the  highest  values  on  the  North- 
European  shelf.  Previous  reconstructions  of  the  surface  currents  in 
the  Iroise  Sea  have  never  documented  velocities  close  to  4  m/s,  nei¬ 
ther  provided  detailed  spatial  patterns  of  velocity  fields.  The  angular 
resolution  of  the  BF  method,  routinely  used  in  processing  the  HFR 
data,  is  too  coarse  and,  in  our  opinion,  unable  to  capture  the  fine 
scale  structure  of  the  surface  flow,  especially  around  islands  and  in 
the  Fromveur  Strait  which  are  of  particular  interest.  An  effort  to  em¬ 
ploy  the  OMA  technique  (Kaplan  and  Lekien,  2007)  for  interpolating 
radial  velocities  provides  much  less  details  (cf.  Fig.  8,  Muller  et  al., 
2009)  due  to  intrinsic  limitations  of  the  method  and  coarse  resolution 
of  radial  velocity  fields. 

Our  analyses  have  shown  that  surface  circulation  in  the  Iroise  Sea 
is  dominated  by  tides  and  winds.  The  wave  induced  component  of 
surface  currents  is  less  important,  while  the  buoyancy  forcing  (due 
to  freshwater  input)  was  negligible  during  the  study  period.  Complex 
bottom  topography,  numerous  islets  and  rocks,  partially  exposed  at 
low  water,  make  the  local  circulation  noticeably  different  from  that 
observed  by  oceanographic  coastal  radars  in  other  regions  of  the 


Fig.  13.  (a)  South-north  component  of  de-tided  surface  residuals  along  the  W-E  line  (shown  by  grey  shading  in  Fig.  12),  and  long-term  time-dependent  surface  residuals  (South- 
north  component)  along  the  W-E  line  (b)  for  April-May  period  2007.  White  vertical  line  indicates  the  location  of  100  m  isobath  on  the  section. 


North-European  tidal  seas  (Carbajal  and  Pohlmann,  2004;  De  Valk, 
1999;  Matthews  et  al.,  1993;  Sentchev  and  Yaremchuk,  2007).  The 
major  differences  are  the  extremely  large  magnitude  of  tidal  currents, 
their  high  spatial  variability  inducing  a  pronounced  eddy  character  of 
near-shore  circulation  resulting  in  a  complex  permanent  structure  of 
the  residual  velocity  field. 

Refined  resolution  allowed  us  to  detect  and  obtain  quantitative  es¬ 
timates  of  the  following  circulation  features: 

1.  A  prominent  dipole  structure  in  the  vorticity  field  characterized  by 
two  oppositely  rotating  eddies,  with  the  maximum  vorticity  of 
—  1.8  10-4s-1  and  1.5  10_4s_1  for  the  northern  and  southern 
eddies,  generated  on  the  leeward  side  of  the  Ushant  Is.  on  flood 


and  ebb  respectively.  Mechanisms  of  vorticity  production  were  de¬ 
scribed  in  detail  by  Pingree  and  Maddock  (1985);  Orbi  and 
Salomon  (1988).  Vorticity  of  the  opposite  sign  is  locally  generated 
by  both  bottom  friction  effects  and  velocity  gradients  with  respect 
to  the  direction  of  tidal  flow.  Advection  by  current  of  the  negative 
vorticity  during  flood,  and  positive  vorticity  during  ebb,  results  in 
appearance  of  two  distinct  zones  of  opposite  vorticity  (dipole 
type  structure)  on  each  side  of  the  island:  negative  vorticity 
north  of  the  Ushant  Island  and  positive  vorticity  south  of  the  is¬ 
land.  Our  experimental  results  are  in  close  agreement  with  these 
theoretical  considerations.  Analysis  of  the  average  divergence 
field  (results  not  shown)  revealed  several  characteristic  structures 
which  seem  to  be  related  to  tidal  generated  eddies.  In  particular,  a 
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Fig.  14.  Same  as  Fig.  1 1  but  for  August-September  2007. 


bi-pole  structure  of  divergence  with  different  polarity  occurs  west 
of  Ushant  Island,  and  is  obviously  related  to  eddies  generated  by 
oscillating  tidal  currents. 

2.  Spatial  intermittence  in  polarization  of  ellipses  has  been  diag¬ 
nosed.  Positive  (ccw)  rotation  of  the  tidal  currents  with  ampli¬ 
tudes  up  to  1.5  m/s  in  the  shallow  water  region  between  the 
Brittany  coast  and  Molene  Islands  has  been  detected.  This  feature 
is  likely  to  be  generated  by  the  bottom  friction  effect,  particularly 
important  at  ebb,  when  a  large  area  between  the  Molene  Islands 
and  the  Brittany  coast  becomes  free  of  water.  This  is  known  as 
“beach  drying”  phenomenon  (Prandle,  1982).  Positive  polarization 
south  of  the  Ushant  Is.  is  related  to  positive  vorticity  (0.5-1. 5 
10“ 4  s_1)  of  currents  in  this  area. 

3.  The  spatial  pattern  of  the  residual  currents  can  be  characterized  by 
two  distinct  zones  approximately  separated  by  the  100  m  isobath: 


In  the  western  (offshore)  zone,  residual  currents  have  a  significant 
contribution  of  the  wind-driven  component,  whereas  the  near¬ 
shore  zone  is  characterized  by  extremely  strong  (up  to  0.4  m/s) 
tidally  generated  residual  circulation. 

4.  Residual  currents  in  the  nearshore  (eastern)  zone  form  two  anticy- 
clonic  eddies:  north  of  the  western  extremity  of  the  Sein  archipel¬ 
ago,  and  north  the  Ushant  Is.  The  maximum  velocity  in  these 
eddies  reaches  0.4  m/s.  The  mean  velocity  field  inferred  from  the 
HFR  data,  contains  a  signature  of  two  along-shore  jets  approxi¬ 
mately  following  the  100  m  (minor  jet)  and  50  m  depth  contours 
(major  jet)  and  merging  south  of  the  Ushant  Is.  in  a  0.5  m/s  strong 
flow.  Qualitatively,  our  results  are  consistent  with  estimates  of  the 
Eulerian  residual  velocities  from  the  numerical  experiments  of 
Muller  et  al.  (2009)  and  from  the  analysis  of  currents  from 
ARGOS  drifters  (Muller  et  al.,  2010). 
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5.  An  extremely  strong  (four-fold)  fortnightly  modulation  of  tidal 
currents  is  documented  northwest  of  the  Ushant  Is.  The  interfer¬ 
ence  of  the  major  semi-diurnal  tidal  constituents  (M2,  S2  and  N2) 
provides  the  major  part  of  the  variability  of  currents  there.  We 
also  diagnozed  a  significant  contribution  of  the  higher  order  non¬ 
linear  tidal  harmonics  (M6,  2MS6,  2SM6)  to  the  surface  current  in 
the  Fromveur  strait.  Here,  the  sixth-diurnal  and  quarter-diurnal 
tidal  harmonics  are  comparable  in  strength  (Table  1),  and  their  in¬ 
terference  affects  the  shape  of  fortnightly  modulation  of  the  cur¬ 
rent  by  flattening  its  maxima  and  amplifying  its  minima. 

A  long-term,  time-dependent,  component  of  residual  currents  has 
been  estimated  and  assessed.  Specific  features  of  this  component  are: 
separation  of  the  basin  in  two  zones,  offshore  and  nearshore,  with  dif¬ 
ferent  control  of  the  RC  and  low-frequency  (weekly)  modulation  of 
the  flow  in  April-May,  2007.  The  comparison  of  the  RC  fields  for 
spring  and  late  summer  periods  revealed  substantial  difference  in 
the  offshore  and  similarity  in  the  near-shore  circulation  patterns  sug¬ 
gesting  a  quasi-permanent  character  of  the  residual  flow  of  the  latter. 

Our  analysis  have  also  shown  that  the  magnitude  of  variation  of 
the  offshore  flow  and  the  mean  velocity  of  permanent  (nearshore) 
flow  are  comparable.  Therefore,  periodic  tidal  forcing  superimposed 
on  the  residual  velocity  shear  may  result  in  strong  variations  of  the 
tracer  trajectories  released  at  different  phases  of  the  tidal  cycle 
(Foreman  et  al.,  1992;  Zimmerman,  1986),  and  cause  strong  disper¬ 
sion  of  the  material. 

The  acquired  data  and  the  presented  results  could  be  very  useful 
for  coastal  circulation  model  validation  and  initialization,  comprehen¬ 
sion  and  interpretation  of  remote  sensing  images  (Le  Boyer  et  al., 
2009),  and  also  for  studying  specific  issues  such  as  eddy,  thermal 
tidal  front  dynamics,  high  resolution  Lagrangian  circulation  and  pas¬ 
sive  tracer  dynamics  in  this  region. 
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Appendix.  Wave-driven  contribution  to  radar  derived  surface  current 

The  surface  current  measured  by  radar  writes: 

U  =  UE  +  USS,  (Al) 

where,  UE  is  the  quasi-Eulerian  surface  current,  i.e.  the  current  at  a 
depth  ofe\/8n  with  \  the  radar  wavelength  (Stewart  and  Joy, 
1974),  and  Uss  represents  the  specific  contribution  of  the  sea  state. 
Let  us  define  a  function  u(f)  as: 

U(f)  =  4n  jg  fo” f  l(E(f ,  0)df  d6,  (A2) 

with /the  wave  frequency,  k  the  wavenumber  and  E(f,Q)  the  sea  sur¬ 
face  elevation  spectrum.  Uss  writes  as  (Broche  et  al.,  1983): 

Uss  =  u(fB)+M ,  (A3) 

where /B  is  the  Bragg  frequency  (0.359  Hz  for  a  radar  frequency  of 
12.4  MHz).  Computations  show  that  M,  which  is  not  developed 


here,  is  a  small  term  compared  to  u  at  present  radar  frequency. 
Then  USs  is  nearly  equal  to  the  contribution  of  the  total  surface  Stokes 
drift,  u(°°),  of  the  waves  with  frequencies  below  fB. 

As  an  output,  the  wave  model  WW3  provides  an  approximation  to 
u(°°),  namely  u(/c),  where  fc  =  0.72  Hz  is  the  cutoff  frequency  of  the 
model.  Recently,  Ardhuin  et  al.  (2009)  proposed  a  robust  approxima¬ 
tion  for  u(f)  (the  norm  of  u(f))  in  terms  of  the  local  10-meter  wind 
speed  Uw  and  the  significant  wave  height  Hs 

u(f)  =  5-10_4(l.25-0.25(0.5//)1-3)l/10min(l/1o,  14.5) 

+  0.025(HS— 0.4)  (A4) 

Eliminating  U10  between  the  two  equations  formed  by  Eq.  (4)  with 
/=/B  and  f=fc  respectively,  and  neglecting  M,  one  can  obtain 

=  0.79  u(fc)  +  5.251 0-3  (Hs— 0.4).  (A5) 

The  analysis  of  the  model  data  for  the  experimental  period  has 
shown  that  (A5)  can  be  further  simplified  to  Lks  =  0.85u(/C). 

We  take  as  direction  of  USs  the  direction  of  u(/c). 
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